Nuclei of Double-Charm Hyperons 
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The ground states of double-charm hyperons form a spin 1/2 isospin 1/2 multiplet analogous to 
that of nucleons. Their main strong interaction may be inferred directly from the corresponding 
nucleon-nucleon interaction by multiplication of the interaction components by the appropriate frac- 
tional difference between interaction strengths for pairs of light flavor quarks and pairs of triplets, e.g. 
nucleons, of light flavor quarks. By construction of the interaction between the recently discovered 
double-charm hyperons by this method from several realistic nucleon-nucleon interaction models it 
is shown that double-charm hyperons are likely to form bound (or possibly meta-stable) states akin 
to the deuteron in the spin triplet state. Double beauty baryons would form corresponding deeply 
bound states. Nucleons and double charm (beauty) hyperons will also form bound states. The 
existence of hypernuclei with double-charm and double-beauty hyperons, which are stable against 
the strong decay, is very likely. 



The SELEX collaboration has recently discovered sev- 
eral double-charm hyperons, the lowest one of which is 
the ground state multiplet with the S++ at 3.46 GeV/c 2 
The ground state multiplet S+ , S+ + forms a spin 
1/2 isospin 1/2 multiplet, with the valence quark config- 
uration dec and ucc |3j. 

The main color-neutral strong interaction between 
double-charm hyperons is due to their light flavor quark 
component. This may be inferred from the nucleon- 
nucleon interaction, by multiplication of the components 
of the nucleon-nucleon interaction by appropriate frac- 
tional coefficients, which relate the interaction strength 
between pairs of light flavor quarks to that between such 
triplets. The interaction between the charm quark pairs 
in different hadrons is expected to be weaker than that 
between light flavor quarks in different hadrons, as the 
latter either arises from the short range interaction that 
is mediated by the exchange of charmonia or from the 
color van der Waals interaction. Because of the isospin 
1/2, the S cc interactions include a long range pion ex- 
change component, in contrast to the A+, which has 
isospin 0. A priori hypernuclei with double charm hy- 
peron should therefore be more likely than hypernuclei 
with single charm hyperons |4|. 

Although the strong interaction between double-charm 
hyperons is weaker than that between nucleons, this is 
partially compensated by their larger mass which weak- 
ens the repulsive effect of their kinetic energy. The latter 
may be schematically illustrated by a calculation of the 
binding energy of the deuteron-like state, which is ob- 
tained by replacement of the nucleon mass by the mass of 
the S cc (3.46 GeV/c 2 ), with the realistic nucleon-nucleon 
interaction models in Refs. 0,0,0, all of which repro- 
duce the experimental binding energy of the deuteron 
and provide a quantitative description of nucleon-nucleon 
scattering data. The corresponding binding energies are 
-71 MeV, -75 MeV and -58 MeV respectively (the scat- 
ter between these values reflects the differences in the de- 
tailed behavior of these interaction models at short range, 



which is not constrained by nucleon-nucleon scattering 
data). 

In operator form the nucleon-nucleon interaction is 
given in terms of rotational invariants of spin and isospin 
as well as momenta and angular momenta. The quark 
model scaling factors between the matrix elements of the 
spin-isospin invariants for S cc and nucleon states may be 
derived from the quark model matrix elements of light 
flavor quark operators ||: 
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Here N represents the nucleon. 

The interaction between two double-charm hyperons 
(and that between double-charm hyperons and nucleons) 
that arises from the interaction between the light flavor 
quarks may be determined from any realistic nucleon- 
nucleon interaction model, that is given in operator form. 
We here consider the models of Refs. 0, 0, mentioned 
above. From these, after dropping the (small) flavor sym- 
metry breaking terms, which are inapplicable to hyper- 
ons (from those which include such), the corresponding 
interactions between double-charm hyperons may be de- 
rived by application of the appropriate downscaling of the 
strengths of the corresponding interaction components. 

In Table 0] the scaling factors for all the NN operators 
that appear in the NN potential models considered are 
listed. Where necessary, the larger scaling factor corre- 
sponding the operator has been taken. The main quali- 
tative difference between the nucleon-nucleon interaction 
and that between double-charm hyperons is the weaken- 
ing by the factor 9 of the central spin and isospin inde- 
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Argonne V18 [5] 
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TABLE I: Quark scaling factors for the components of the 
interaction models of the form 0, [1| . 



pendent component, which in the case of nucleons pro- 
vides both the intermediate attractive interaction and the 
strong short range repulsion. The fact that the strengths 
at short range of the components of the three nucleon- 
nucleon interaction models differ considerably is of lit- 
tle numerical significance for the description of nucleon- 
nucleon scattering observables at low energy, because of 
the strong short range repulsion. The weakening of that 
repulsive interaction component in the double-charm hy- 
peron interaction increases the significance of the differ- 
ences in the short range interaction components, which 
is revealed as large differences in the predicted binding 
energies. 

For the deuteron the main cause of binding is the 
strong isospin dependent tensor force, the long range part 
of which is due to pion exchange. Since this is weaker by 
a factor 25 for double-charm hyperons, it follows that 
its role is significantly smaller for double charm hyper- 
ons. Below we show, on the basis of the quark model 
scaling factors, that two of the three rescaled nucleon- 
nucleon interaction models nevertheless do imply that 
double-charm hyperons form deuteron-like bound states, 
with binding energies in the range 87 — 457 MeV. 

The binding energies were calculated by diagonaliza- 
tion of the Hamiltonian in a discretized Gauss-Legendre 
mesh in momentum space. The bound state equation can 
be written as, 
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Here I and I' are orbital angular momenta, [i is the re- 
duced mass of the system and Vu'(q,q') is the partial 



wave decomposed potential (< qlsj \V\q'l'sj >). For 
the purpose of the present investigation a mesh of 80 
Gauss-Legendre points was sufficient for stable results. 
As a test of the numerical method the deuteron binding, 
—2.2 MeV, was calculated with all the models, AV18, 
Nijmegen II and Paris potential. 

In Table [H] the calculated binding energies obtained 
for the deuteron-like states of double-charm hyperons are 
given. The scatter between the calculated values provides 
an estimate of the theoretical uncertainty that derives 
from the different short range behavior of the nucleon- 
nucleon interaction models. 

It is natural to assume that the experimentally dis- 
covered double-charm hyperon state S+ (3520) represents 
the double-charm spin 3/2 state, which is analogous to 
the A(1232) [103- The 60 MeV splitting between this 
state and the S++ (3460) is in line with quark model esti- 
mates ^]], although slightly smaller than numerical lat- 
tice method based estimates As the coupling to 
this state can only increase the calculated binding en- 
ergy, the binding energies in Table [I] may be viewed as 
lower bounds. 

Two-baryon states formed of double-charm hyperons 
can in principle couple to states with a single charm and 
a triple-charm Slccc by the quark rearrangement inter- 
action. If the latter states have lower energy the former 
are metastable rather than bound. This depends on the 
size of the binding energy as compared to the mass dif- 
ference: 
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Here I represents a light quark. For many non-relativistic 
quark models the inequality A c < holds 0. Corre- 
sponding numerical estimates suggest that A w [130 — 
158] MeV 01- Adoption of those values imply that the 
double charm hyperons form bound states with the AV18 
potential, but only metastable states with the Nijm II 
potential. 

To further explore the theoretical uncertainty of the 
calculated binding energies these have also been calcu- 
lated using the class "AVn" of systematically simplified 
versions of the V18 5] interaction models [J. These re- 
sults are shown in Table ITTTl The number of each poten- 
tial designates the number of operator structures present, 
ordered as in the first column of Table [i] These poten- 
tial models are reprojections of the full AV18 potential, 
of which all but AVI' and AVX' reproduce the deuteron 
binding energy. 

The results show a considerable spread in the calcu- 
lated binding energies. No bound state is found with 
the AV8' potential, which already includes tensor and 
spin orbit forces. With the simpler AV2', AVX', AV4' 
and AV6' models bound states appear. Consequently 
nucleon-nucleon interaction models, which give the cor- 
rect deuteron binding energy, with (AV6') and without 
tensor or spin-orbit forces, predict a bound state in the 
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Double-Charm hyperons 


Potential 


Binding Energy (MeV) 


AV18 


-457 (-28) 


Paris 




Nijm II 


-87 


Double-Beauty double-charm hyperons 


AV18 


-603 (-183) 


Paris 


-0 


Nijm II 


-102 


Double-Beauty hyperons 


AV18 


-782 (-439) 


Paris 


-2 


Nijm II 


-123 (-20) 



TABLE II: Binding energies for the 5i + — S^ c and 5® b — E hb 
systems obtained with Argonne V18 |3, AV18, Nijmegen II |6j 
and Paris potentials. For the double-beauty baryons the 
mass values: Mbbu — Mbbd ~ 9000 MeV were used. The value 
in brackets corresponds to a second bound state. 
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deuteron 


Potential 


Binding Energies (MeV) 


AV8' 




-0 


-3 




-2.2 


AV6' 


-1 


-5 


-15 




-2.2 


AV4' 


-30 


-41 


-58 (-4) 




-2.2 


AVX' 
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-22 


-33 




-0.4 


AV2' 


-26 


-48 


-88 (-4) 




-2.2 


AVI' 










-0.4 



TABLE III: Estimated binding energies for the pairs H cc — E cc , 
H cc — Hfci, and Ef,i, — of double-heavy together with the 
deuteron as obtained with the class of Argonne potentials in 
Ref. ®. 

S+ + — 5+ system in the spin triplet state. The case of 
the AV2' potential, for which the quark model scaling is 
simply a factor 1/9 of the 1 part, and which predicts a 
bound state of 26 MeV, is particularly telling. The bind- 
ing energies obtained with the AVn' class of potentials 
are considerably smaller than the one obtained with the 
full AV18 potential. In particular, with the numerical 
estimates in Ref [l3j ] based on the mass inequalities, see 
Eq. 0, these bound states would be metastable against 
decay to the single-charm - triple-charm hyperon state. 

The origin of the large binding energy given by the 
AV18 interaction model may be traced to its large isospin 
dependent squared spin-orbit interaction, which acts in 
the D— state. The relative significance of this interac- 
tion component derives from the large quark model scal- 
ing factor (Table 1). The general reason for the large 
spread in the calculated binding energies obtained with 
the different interaction models is the fact that their 



components differ notably at ranges less than 0.6 fm. 
These differences are of little significance for nucleon- 
nucleon scattering observables at low energy, because of 
the very strong short range repulsion. The reduction of 
the strength of the repulsive central interaction by a fac- 
tor 9 enhances the role of the short range behavior of 
the state dependent components of the interaction in the 
case of double-charm hyperons. 

Double beauty-hyperons analogous to, but heavier 
than, the double-charm hyperons are expected to exist in 
the mass range above twice the B-meson mass. The in- 
teraction between the ground state multiplets of double- 
beauty hyperons with the quark configurations (bbu) and 
(bbd) is expected to be very similar to that of the double 
charm hyperons, as their main interaction is that between 
their light flavor quarks. Their binding energy will how- 
ever be much larger than that of double-charm hyperons 
in view of their larger mass. This is illustrated in Table UTI 
where the binding energy of two double-beauty hyperons 
is estimated by the same method as used above for dou- 
ble charm hyperons with the assumption that their mass 
is: M bbu = Mbbd « 2m fc « 2 x 4500 MeV. In Table |ni| 
the binding energies obtained with the AVn' potentials 
are listed for comparison as well. The majority of these 
potential models yield a bound state, in some cases two, 
with binding energies up to few hundred MeV. 

The possibility of metastability is in this case bounded 
by a similar inequality, as in the case of the charmed 
hyperon case: A b = Mbbb + Mui — 2Mbu < 0, where 
its estimated value ranges [348 - 372] MeV [H Q As 
in the case of the double-charm hyperons the result ob- 
tained with the AV18 potential is a bound state, while 
that obtained with the Nijm II potential is a metastable 
state. 

Finally we have also explored the possibility of 
deuteron-like bound states of nucleons and double-heavy 
hyperons: N — E cc and N — Ebb- Such bound states were 
found with the AV18 and the Nijm II potential models. 
The AV18 potential gives bound states at —388 MeV and 
—494 MeV for the N — E cc and N — Ebb systems respec- 
tively. The Nijm II potential gives bound states at —35 
MeV and —76 MeV for the N — E cc and N — Ebb systems 
respectively. 

In summary it has been shown that bound states of the 
recently discovered double-charmed hyperons are likely. 
Their binding energies were estimated by construction of 
their main strong interaction by the quark-model scaling 
relations from several realistic phase-equivalent nucleon- 
nucleon interaction models. The E cc hyperons are likely 
to be bound in the spin triplet state by more than 30 
MeV. The Ebb hyperons are expected to be bound by 
more than 50 MeV in the triplet state, while pairs of E cc 
and Ebb baryons are expected to be bound with binding 
energies of that order. The qualitative nature of these 
estimates reflects the uncertain short range behavior of 
the extant realistic nucleon-nucleon interaction models. 
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